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Two strongly-pumped parametric interactions are simultaneously realized in a single nonlinear 
crystal in order to generate three strongly correlated optical fields. By combining together the 
outputs of two of the three detectors measuring intensities of the generated fields, we obtain the 
joint photocount statistics between the single field and the sum of the other two. Moreover, we 
develop a microscopic quantum theory to determine the joint photon-number distribution and the 
joint quasi-distributions of integrated intensities and prove nonclassical nature of the three-mode 
state. Finally, by performing a conditional measurement on the single field, we obtain a state 
endowed with a sub-Poissonian statistics, as testified by the analysis of the conditional Fano factor. 
The role of quantum detection efficiencies in this conditional state-preparation method is discussed 
in detail. 
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I. INTRODUCTION 

The quantum properties of photon pairs generated in 
theprocess of spontaneous parametric down-conversion 
[H 09 have been investigated in many experimental and 
theoretical works during the last 30 years. Pairwise char- 
acter of the fields generated in this process has been used, 
ej*., for testing fundamental laws of quantum mechanics 
[If, as well as for quantum teleportation Q, quantum 
cryptography || and in metrology applications [6] . The 
theory that describes this kind of interaction has been 
elaborated from several points of view for fields contain- 
ing just a fraction of a photon pair 0, H, Q and for fields 
composed of many photon pairs [10| . Also stimulated 
emission of photon pairs has been addressed [ll|, EH EH • 

More recently, the development of the field of 
quantum-information processing [14j | has drawn attention 
to three-field quantum correlations. In fact, interesting 
correlations can be reached when strongly-pumped para- 
metric down-conversion and parametric amplification are 
combined together through a common field. Such a sys- 
tem can be built in a single nonlinear crystal oriented 
in such a way that phase-matching conditions for both 
the interactions are fulfilled together [l5| . It has been 
demonstrated that the tripartite state generated in this 
way is endowed with entanglement in the number of pho- 
tons. In particular as the constant of motion admitted 
by the hamiltonian that describes the process suggests, 
we obtain that the number of photons in one of the three 
generated fields is always equal to the sum of the pho- 
tons in the other two fields These properties make 
the system useful for the generation of nonclassical states 
by means of a suitable conditional state preparation. In 
more detail, if a given number no of photons in field ao 
is detected then the remaining two fields a\ and 02 are 
ideally left in a state \ip) = Ya=o c zlOil n o — O2 entangled 



in photon numbers, \l)i means the Fock state with I pho- 
tons in field and coefficients q depend on the nonlinear 
interaction. 

Here we present the experimental realization of this 
scheme and in particular we study the joint photocount 
statistics between one field (single field) and the sum of 
the other two (compound field) in order to find out the 
experimental conditions in which it is possible to realize 
a deterministic source of states entangled in the number 
of photons. We note that this source re quir es photon- 
number resolving detectors 0, El, El HI 1M HI . 

Note that fields composed of photon pairs have been 
experimentally investigated [UHliHs, 26, 27] under con- 
ditions that allowed having up to several thousands of 
photon pairs per pump pulse. A detailed theory based 
on a multi-mode description of the generated fields has 
been developed for spontaneous [28j as well as for stim- 
ulated processes [29[ with the aim to interpret the ex- 
perimental data. In this work, as the photocounts from 
two of the three detectors are combined together, we also 
use this theory in order to correctly interpret and process 
the measurements. In particular, we can determine the 
joint single-compound-field (JSCF) photon-number dis- 
tribution and the JSCF quasi-distributions of integrated 
intensities from the measured JSCF photocount distribu- 
tion. Moreover, the conditional photocount and photon- 
number distributions, together with the corresponding 
Fano factors, are derived in order to point out the versa- 
tility of the scheme as a source of nonclassical states. To 
this aim, the quantum detection efficiency is an impor- 
tant parameter. 

Theoretical description of the nonlinear process is pre- 
sented in Sec. II. Photon-number distributions and quasi- 
distributions of integrated intensities as well as other 
properties of the measured fields are derived in Sec. III. 
Experimental results are discussed in Sec. IV. Sec. V 
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gives conclusions. 

II. THREE-MODE PARAMETRIC PROCESS 

We consider a three-mode parametric process in which 
two strong coherent classical fields pump two interlinked 
nonlinear interactions: a frequency down-conversion and 
a parametric amplification. The correspondin g in terac- 
tion Hamiltonian can be written as follows (TBI l30j] : 

H^t = 7oOo a l + lia\a2 + h.c, (1) 

where a, i = 0,1,2, are the corresponding photon 

annihilation (creation) operators and 70 and 71 are cou- 
pling constants for frequency down-conversion and para- 
metric amplification, respectively. Symbol h.c. stands 
for the hermitian conjugated term. Analytical solutions 
of the corresponding Heisenberg equations for the annihi- 
lation and creation operators are given in [1 51 ] . They can 
be used for the determination of the normal three-mode 
characteristic function of the spontaneous process: 

C(fa,fa,fa) = exp[-A)|/?o| 2 -A|/3i| 2 - 5 2 |/?2| 2 
+(D 01 l%ft + D m ft% + D 12 fa[3* 2 + c.c.)], (2) 

where /3j (i = 0,1,2) are parameters, c.c. means the 
complex conjugated term and 



Bo 


= (Aa^Aao) = 


l/i| 2 + l/ 2 


A 


= (AajAai) = 


l5o| 2 , 


B 2 


= (Aa^Aa 2 ) = 




An 


= (AaoAai) = 


fo*9o, 


An 


= (Aa Aa 2 ) = 


foho, 


A2 


= -(Aa\Aa 2 ) 


= -ho9o- 



The functions go, ho, and /o are defined in [la ]. In order 
to study the nonclassical nature of the three-mode state, 
we introduce the determinants 



Ki3 


= B\B 2 — 


\D 12 \ 2 =0, 






K01 


= BqBi — 


|A»i| 2 = - 


\9o\ 2 


<o, 


K 02 


= b b 2 — 


|Ab| 2 = - 


\ho\ 2 


< 0. 



According to Eqs. ([JJ fields a,\ and a 2 are classically cor- 
related whereas correlations between fields ao and a% (clq 
and 02) can lead to nonclassical behavior. 

In the experiment, the outputs of the detectors placed 
on fields a\ and a 2 have been summed with the aim to 
measure photocount correlations between the single field 
ao and the compound field formed by fields a\ and a 2 . 
These correlations are important to test the performance 
of a source of states entangled in the number of photons 
and obtained after a conditional measurement of no pho- 
tons in field ao- They can be derived from the following 
normal characteristic function: 

C(fa, fa, fa) = exp[-B |/3 | 2 - B 12 \fa\ 2 

+ (D 0A2 p* Pl+h.c.)l (5) 



where 

B 12 = B x +B 2 - D 12 - D 21 = \ho + .go 1 2 , 

A),i2 = An + Aj2 = / *.9o + foho, 

Ko,i2 = BoB l2 -\DoM 2 = -\go\ 2 -\ho\ 2 < 0. (6) 

To obtain Eq. ([S]) we have assumed that -fi is real. We 
note that information about the coupling constants 70 
and 71 can be obtained from the reconstruction of the 
photocount distributions provided that the interaction 
time t is known. 

We assume that each field is composed of M indepen- 
dent temporal modes [27j • Then the variances of inte- 
grated intensities W can be expressed as follows: 

((AW0) 2 ) = MB 2 , 

((AW,) 2 ) = MB],j = 1,2, 

(AW AWj) = M\D 0j \ 2 ,j = 1,2, 

(AW!AW 2 ) = M\D 12 \ 2 . (7) 

If photocounts belonging to fields a\ and a 2 are combined 
together, we can write the variances of the integrated 
intensities W in the following form: 

((AW 12 ) 2 ) = ((AW 1 ) 2 ) + ({AW 2 ) 2 )+2(AW 1 AW 2 ) 

= M\B 12 \ 2 = M(B 2 + B 2 + 2\D 12 \ 2 ) 

= A1(B 1 +B 2 ) 2 , 

(AW AW 12 ) - Af|A,i2| 2 = M|/o| 2 (lffo| 2 + l^o| 2 )- (8) 

Real value of 71 has been again assumed. 

Note that the above-presented analysis can be general- 
ized to stimulated processes [3l|. However, single as well 
as compound fields have to be stimulated in order to sup- 
port nonclassical effects by interference terms. As for the 
compound field, stimulation of one of its components (i.e. 
field ai or a 2 ) is sufficient to observe increased nonclas- 
sical effects. Stimulation of only one field (i.e. field a Q or 
a± or a 2 ) results in increased values of noise only. 

III. PHOTON-NUMBER DISTRIBUTIONS AND 
QUASI-DISTRIBUTIONS OF INTEGRATED 
INTENSITIES 

The experimental scheme used to generate two non- 
linear interactions is depicted in Fig. [T] The harmonics 
of a continuous-wave mode-locked Nd:YLF laser regen- 
eratively amplified at a repetition rate of 500 Hz (High 
Q Laser Production, Hohenems, Austria) provided two 
pump fields. In particular, the third harmonic pulse at 
349 nm (~ 4.45 ps pulse-duration) was exploited as the 
pump field a p o in frequency down-conversion, whereas 
the fundamental pulse at 1047 nm (~ 7.7 ps pulse- 
duration) was used as the pump field a p \ in parametric 
amplification. The two interactions simultaneously satis- 
fied energy-matching (lu p o =^0 + ^2 an( ^ ^1 = ^2 + ^pi ) 
and type I phase-matching (kSj = kg+kSJ, kj = kj+kM 
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FIG. 1: Scheme of the experimental setup: BBO, non- 
linear crystal; NF, variable neutral-density filter; Pi, pin- 
holes and Di, p-i-n photodiodes, i — 0,1,2; fi, lenses, 
i — 0, 1, 2,pQ,pl, pi'; PRE+AMP, low-noise charge-sensitive 
pre-amplifiers followed by amplifiers; SGI, dual-channel syn- 
chronous gated-integrator; ADC+PC, computer integrated 
digitizer. 



conditions, in which ujj are the angular frequencies, k, 
denote the wave vectors and suffixes o and e indicate 
ordinary and extraordinary field polarizations. As de- 
picted in the inset of Fig. [TJ we set the pump-field a p o 
direction so that the wave vector k p o was normal to the 
crystal entrance face and propagated along the z-axis of 
the medium. We also aligned the wave vector k p i of the 
other pump field a p \ in the plane {y, z) containing the 
optical axis (OA) of the crystal and the wave vector k p o • 
As the nonlinear medium we used a /3-BaB2C>4 crystal 
(BBO, Fujian Castech Crystals, China, 10 mm x 10 mm 
cross section, 4 mm thickness) cut for type-I interaction 
(^cut = 38.4 deg), into which both pumps were strongly 
focused. Typical intensity values of the pump fields were 
~ 5 GW/cm 2 for a p0 and ~ 2 GW/cm 2 for a p i. The re- 
quired superposition in time of the two pump fields was 
obtained by a variable delay line. 

As we have already shown in Ref. we decided 

to generate three fields at non-degenerate frequencies 
by choosing a phase-matching condition in the plane 
(y, z) [32]. In order to investigate the nature of the 
state obtained by the interlinked interactions, we se- 
lected a triplet of coherence areas by means of pin-holes, 
whose sizes and distances from the crystal were chosen 
by searching for the condition of maximum intensity cor- 
relations between the generated fields [33j. In fact, we 
expect strong correlations not only between the number 
of photons in the field ao and the sum of the other two 
fields (compound field), but also singularly among the 
numbers of photons in all pairs of fields. By applying 
this criterion, we put two pin-holes of 30 fiia diameter 
at distances do = 60 cm and di = 49 cm from the BBO 
along the path of the signal beam at 632.8 nm and of the 
idler beam at 778.2 nm, respectively. Moreover, as the 



beam at 446.4 nm has smaller divergence compared to the 
other two fields, we selected it by means of a 50 /an diam- 
eter pin-hole placed at a distance d\ = 141.5 cm from the 
crystal. The light was suitably filtered by means of band- 
pass filters and focused on each detector. In particular, as 
we performed measurements in the macroscopic intensity 
regime (more than 1000 photons per coherence area), we 
used three p-i-n photodiodes (two, Do and Di in Fig. [TJ 
S5973-02 and one, D 2 , S3883, Hamamatsu, Japan) as the 
detectors. We obtained the same overall detection effi- 
ciency (rj — 0.28) on the three arms by inserting two ad- 
justable neutral-density filters in the pathways of fields a\ 
and a2. The current output of the detectors was amplified 
by means of two low-noise charge-sensitive pre-amplifiers 
(CR-110, Cremat, Watertown, MA) followed by two am- 
plifiers (CR-200-4 us, Cremat): to this aim, we connected 
the detectors Di and D 2 to the same amplifier device by 
means of a T-adapter. The two amplified outputs were 
then integrated by synchronous gated-integrators (SGI 
in Fig. [TJ SR250, Stanford Research Systems, Palo Alto, 
CA) sampled, digitized by a 12-bit converter (AT-MIO- 
16E-1, DAQ National Instruments) and recorded by a 
computer. 

From the collected experimental data we obtained the 
first ((mo), (toi 2 )) and the second ((m 2 ,), {m\ 2 )) mo- 
ments of the photocount (photoelectron) distribution of 
the single and compound fields, respectively. Moreover, 
the correlation between the number of photoelectrons in 
the single field and that in the compound field was mea- 
sured. The unavoidable contributions of additive noise 
present in the detection chain can be quantified in an 
independent measurement and then subtracted from the 
experimental data. By correcting the moments of the 
photoelectron distribution for the quantum detection ef- 
ficiency, we can obtain the moments for photons. Here 
we present the first ((n ), (nu)) and the second ( (no), 
(n 2 2 ), (non^)) moments of the photon- number distribu- 
tion: 



(rij) 
(n 2 > 
(n ni 2 ) 



(TO 2 )/77 2 -(l-77)K)/n 2 , z = 0,12, 



(TO mi 2 )/?7 2 . 



(9) 



In addition, the moments for the integrated intensities 
can be determined as follows: 



(Wi) = (n,), 

(IV, 2 ) = (n 2 )-(n 4 ),z = 0,12, 
(WqWu) = (n ni 2 ). 



(10) 



In a theory based on the generalized superposition of 
signal and noise properties of the fields can be quantified 
using coefficients Bq, Bi 2 , -Do, 12 and number M of tem- 
poral modes. These quantities can be determined from 
the moments of integrated intensities [28| : 



B, 

Mi 

|A),12| 



((AW^/iWi), 
(Wrf/dAWt) 2 ), i = 0,12, 
^(AW AW 12 )/M. 



(11) 
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The mean number of photons in mode i is given by Bi 
(i = 0, 12) whereas D 12 quantifies the correlations be- 
tween the single field and the compound one. The num- 
ber of modes M can be determined cither from the exper- 
imental data measured in the single or in the compound 
field (see Eqs. (fTTj)). Ideally, Mo should be equal to M12 
(Til l . This cannot be reached in a real experiment because 
of non-perfect alignment and detection noise. However, 
a correct alignment of the experimental setup allows to 
have M « M 12 so that we can define M = (M +M 12 )/2. 
A more detailed analysis concerning the determination of 
M can be found in [34| . 

The JSCF photon-number distribution p(no,n\ 2 ), 
which can be derived from the normal characteristic func- 
tion C in Eq. ([5]), is written as follows (H)]: 

, v 1 (Bo + K . 12 ) n »{B 12 + K , 12 ) n ^ 

P{n a , U 12 ) r(M) (1+Bq + B ^ + Ko l2) n 0+ n 12 +M 

min(n ,ni2) r , , . >, n 

v T(n a + n l2 + M — rj 

^ r\(n - r)\(ni2 - r)\ 

(-Ko,i2) r (l + B a + B 12 + K , 12 ) r 
X [(Bo + K .i 2 )(B 12 +K OA2 )Y ' [ > 

where T is the gamma function. If determinant K 12 
defined in Eqs. ((6]) is negative, the given field cannot 
be described classically. The quantities Bo + Kq.\2 and 
B12 + if 0,12 occurring in Eq. (fT2| cannot be negative 
and can be interpreted as components of fictitious noise. 
In the ideal lossless case, #0,12 = ~Bq = —B\ 2 and then 
the joint photon-number distribution p(no, rc.12) takes the 
form of the diagonal Mandcl-Rice distribution. As fre- 
quency down-conversion produces couples of photons, 
we expect that the highest values of the elements of 
p(no,ni 2 ) are near the diagonal no = n\ 2 . In addition 
certain classical inequalities can be violated in this region 
0. 

If ifo. 12 = 0, i.e. at the boundary between the classical 
and nonclassical behaviors, the compound Mandcl-Rice 
formula for JSCF photon-number distribution p(no, n\ 2 ) 
can be simplified as follows 



p(n ,ni 



T(n +n 12 + M)BZ°B™ 
T(M)n \ nu! (1 + B + B 12 ) n °+ n ™+ M ■ 



(13) 

In the ideal case only the diagonal elements of the 
distribution p(no,n\2) are different from zero: in this 
case the post-selection scheme for the preparation of a 
state entangled in the number of photons perfectly works. 
However, losses present in any experimental implemen- 
tation cause discrepancy from this ideal situation. In 
order to quantify this discrepancy, we determine the 
conditional compound-field photon-number distribution 
Pc,i2{n\2',no) provided that no photons are detected in 
the single field: 



Pc,i2(ni 2 ;n ) 



p(n ,n 12 )/^2p(n ,k). 

k=0 



(14) 



Fano factor -F C ,12 of photon- number distribution 
Pc,i2(ni2',no) can be expressed as follows: 

^,12(^0) = 1 

(1 + M/n )[(B 12 + K 0yl2 )/(1 + Bp)} 2 - {K , 12 /B ) 2 
(1 + M/n ){B 12 + K Q , 12 )/{1 + B ) - K0.12/B0 
w 1 + Jfo.ia/So. (15) 

Note that the last approximation in Eq. (fT5|) holds 
for i^o,i2 ~ —B12 and highlights that negative values 
of determinant -Ko,i2 cause sub-Poissonian conditional 
photon-number distribution. For the ideal lossless case, 
i£o,i2 = —Bo = —B\ 2 and Fano factor F c ^\ 2 is equal to 
0. ' 

The correlations in the number of photons can also be 
quantified by the distribution p_ of the difference photon 
number no — n\ — n 2 defined as: 

00 

P~( n )= X! < 5n,« -n 1 -n 2 P(^0,"-l,"-2), (16) 
no, ni,n2=0 

where 5 is the Kronecker symbol. The variance of the 
difference no — n\ — n 2 can be lower than the sum 
(no + n± + n 2 ) of the mean photon numbers in all fields. 
This happens for nonclassical fields and we obtain sub- 
shot-noise correlations in this case 0, [26j. Under this 
condition the so-called noise reduction factor 1 26] , 



R = 



{[A(n - ni -n 2 )] 2 ) 
(n ) + (ni) + (n 2 ) 



(17) 



is lower than 1. 

The JSCF photon-number distribution p(rtp, ni 2 ) is 
given by Mandel's photo-detection formula [3, [2^, l3a |. 
which is defined in terms of the s-ordered quasi- 
distribution P s (Wo, W\ 2 ) of integrated intensities. This 
formula can be inverted in such a way that the quasi- 
distributions of integrated intensities can be written in 
terms of the distribution p(no,n\ 2 ) in Eq. (|12[) with 
the following results. If we have the s-ordered deter- 
minant iifo,i2« > (#o,i2s = B 0s Bi 2s - \D Q , 12 \ 2 , B iiS = 
Bi + (1 - s)/2, i = 0,12), the s-ordered JSCF quasi- 
distribution P s (Wo,Wi 2 ) of integrated intensities is a 
non-negative ordinary function [281 ] : 



Ps(W ,W 12 ) = 



1 



T(M)K™ 



K 2 12s W W 12 
|A)- 12 



(M-l)/2 



x exp 



12s \ 1-^0,121 

(B 12s Wo/Bq s + W 12 )Bq s 

Kq,12s 



where Im is the modified Bessel function. 

On the other hand, for K 0: i2 S < 0, the JSCF quasi- 
distribution P s (Wo,Wi 2 ) of integrated intensities takes 
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the form of a generalized function. It can be approxi- 
mated using the following expression (28j : 



A(W Wi2) (M - 1)/2 



7rr(M)(i?o s i?i2 S ) M / 2 
W W u 




x sine 



(19) 



in which sinc(x) = sin(x)/x and A = (— .Ko,i2s) _1 • The 
quasi-distribution derived in Eq. (jTHJ) typically oscillates 
and has negative values in some regions. The threshold 
value Sth giving the boundary between the expressions in 
Eqs. (fT5)) and (|19p can be determined from the condition 
K ,i2s = 0: 



Sth 



1 + B + B u 



{Bo + B 12 f 



(20) 



in which — 1 < s t h < 1- 

The variances of the difference n — n 12 of single- and 
compound-field photon numbers (n\i = n± +712) arid the 
difference Wo — W12 of the single- and compound-field 
integrated intensities (W12 = W\ + W2) are linked by the 
following formula: 

([A(n -n 12 )] 2 } = (no) + (n 12 } + (lA(Wo~W 12 )] 2 ). (21) 

Equation (f2"Tj) shows that negative values of the quasi- 
distribution P s (Wo , W12) of integrated intensities are 
needed to observe sub-shot-noise correlations in single- 
and compound-field photon numbers, i.e. R < 1 in 
Eq. dm). 



IV. EXPERIMENTAL RESULTS AND THEIR 
INTERPRETATION 

From the intensity measurements of the single and 
compound fields we calculated the following moments 
for photoelectrons: (m > = 1225.183, (m 12 ) = 1186.138, 
{ml) = 1609827, (m\ 2 ) = 1518257, (m m 12 ) = 1562402. 
Moreover, the knowledge of the quantum detection effi- 
ciencies 770 = rj\ = f]2 — 0.28 allowed us to derive the 
corresponding moments for photons: (no) = 4375.654, 
(ni 2 ) = 4236209, (n 2 ) = 20522260, (rc? 2 ) = 19354630, 
(n n 12 ) = 19928600. 

In the case of photoelectrons the calculated values of 
the coefficients in Eqs. (TTTj) are: 

Bo = 87.765104, £12 = 92.861384, 

I -D0.12I = 90.371968, M = 13.3665. (22) 

We also note that the number of modes Mq = 13.95980 
and M12 — 12.77321 relative to the single and compound 
fields, respectively, are almost the same. In the case of 



photons the values of coefficients Bo, Bi 2 and | -Do, 12 1 are 
the following: 

Bo = 313.447,^12 = 331.648, 

|A),i2 1 = 322.757. (23) 

The number of modes M is the same for photocount and 
photon-number distributions. 

Determinant -Ko,i2 is negative both for photoelectron- 
(i^o,i2 = —17.104) and photon- number (K0.12 = 
—218.158) distributions. At the same time, the fluctu- 
ations in the difference between the single and the com- 
pound fields both for photoelectrons (R = 0.954) and for 
photons (R — 0.837) are below the shot-noise level. In 
accordance with Eq. (|21[) . this means that the wave vari- 
ance ([A(Wo- W 12 )} 2 } of the difference between the single 
and compound field integrated intensities is negative: in 
fact we have the wave variance equal to —110.572 for pho- 
toelectrons and —1406.699 for photons. These negative 
values are due to the fact that frequency down-conversion 
emits the same number of photons into the single field 
do and in the compound field formed by fields a\ and 
a.2- For the same reason we obtain that the covariance 
C = (AmoATOi 2 )/ v / ((Am ) 2 )((Ami2) 2 ) for photoelec- 
trons and the analogous expression for photons assume 
values very close to 1: 0.991 for photoelectrons and 0.990 
for photons. Moreover, we note that even the principal 
squeezing parameter A = 1 + Bo + B12 — 2|£>o,i2 1 
indicates a nonclassical behavior: we obtained A = 0.882 
for photoelectrons and 0.581 for photons; A = 1 holds 
for coherent states. On the basis of the above mentioned 
quantities, we want to emphasize that the quantum na- 
ture of the state produced by the nonlinear process is 
more evident in terms of photons than in terms of pho- 
toelectrons. 

The JSCF photocount, p(mo,mi2), and photon- 
number, p(no,ni2), distributions calculated along the 
formula in Eq. (| 12(1 with values of parameters appear- 
ing in Eqs. flit]) and flUJ are shown in Fig. [2] and Fig. El 
respectively. There are strong photon-number correla- 
tions between the single field and the compound one. In 
particular, we note that the elements of p(mo, TO12) and 
p(no, 7112) assume the highest values in the vicinity of the 
diagonal, i.e. for mo ~ tu\2 and no ~ H12 (see contour 
plots in Figs. and [3]). 

Values of Fano factors F Cj i2 of the conditional pho- 
tocount and photon-number distributions are impor- 
tant from the point of view of the conditional state- 
preparation scheme. We remind that a nonclassical state 
requires < 1, i.e. sub-Poissonian statistics of the 

conditional distributions. In our case this requirement is 
not fulfilled by the photocount distribution, as shown in 
Fig. [H This behavior is due to the relatively low quan- 
tum detection efficiency (77 = 0.28). On the other hand, if 
the experimental data are corrected for the quantum ef- 
ficiency, we can obtain a sub-Poissonian photon-number 
distribution for no > 5 (see Fig. [5]). The greater the 
number n of detected single-field photons the smaller 
the value of Fano factor F c 12. Actually, the values of 
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FIG. 2: JSCF photocount distribution p(mo, m\2) (a) and its 
contour plot (b) 



Fano factor F c ,i2 do not change for numbers of mo and 
no greater than 100. 

The behavior of the JSCF quasi-distribution 
Ps(Wq,Wi2) of integrated intensities depends on 
the value of the ordering parameter s. In our case the 
threshold value of the ordering parameter, sth, is equal 
to 0.811 for photoelectrons. This means that the non- 
classical behavior of the quasi-distribution P a {Wo,Wi2) 
is expected only for values of s greater than sth- Indeed, 
in Fig. [6] we show that oscillations and negative values of 
P s (Wq, W12) occur for s = 0.9 {K 0i i 2s < 0); on the other 
hand, the nonclassical features do not appear for s = 0.4 
(-Ko,i2s > 0). This means that the quantum noise present 
in the detection chain covers the nonclassical behavior. 
In the case of photons, the JSCF quasi-distribution 
P s (Wq, W12) has a lower threshold value of the ordering 
parameter, namely Sth = 0.324. For this reason, the 
JSCF quasi-distribution P s (Wo,Wi2) shown in Fig. [7] 
already shows nonclassical features like oscillations and 
negative values for s — 0.4. 

It is clear from all these considerations that to ap- 
preciate the quantum features of the three-mode state 
and its usefulness for the preparation of a conditional 
state, the quantum detection efficiencies must be high 



FIG. 3: JSCF photon-number distribution p(no, nil) (a) and 
its contour plot (b). 
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FIG. 4: Fano factor F C) i2 of the compound-field conditional 
photocount distribution p c ,i2 as a function of the number mo 
of detected single-field photoelectrons. 



enough. A typical dependence of the conditional Fano 
factor ^,12(^0) on the quantum detection efficiency 77 
(v = ? /o = ?7i = 772) is plotted in Fig. [8] a) for no = 
3000. By decreasing the value of 77 we obtain an in- 
crease in the value of F c _i 2 (no). Moreover, it can be 
demonstrated that there is a threshold value of the effi- 
ciency 7y cr it below which the conditional compound-field 
photon-number distribution is no more sub-Poissonian. 
As shown in Fig. [5]b), the critical value 7y cr it decreases 
by increasing the number no of the photons in the single 
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FIG. 5: Fano factor F Ct \2 of the compound-field conditional 
photon- number distribution p c ,i2 as a function of the number 
no of detected single-field photons. 

a) 
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FIG. 6: JSCF quasi-distributions P s (W ,Wi 2 ) of the single- 
field (Wo) and compound-field (W12) integrated intensities 
corresponding to photoelectrons for s = 0.4 (a) and s — 0.9 
(b). 

field. We note that graphs in Figs. [S] and [9] have been ob- 
tained using the formula for Fano factor F C) i2 in Eq. (fl"5]) 
assuming substitution Bq — > 770-B0, B12 —> t\\Byi, and 
Do 12 ~~ * VoViDq yi and coefficients Bq, B12, and £>o,i2 
for'photons [Eqs.' ([25])] 

In more detail, only the value of the quantum detec- 
tion efficiency t]q is crucial for the sub-Poissonian be- 
havior of the conditional compound-field photon-number 
distribution, as shown in Fig. [9l where the contour plot 
of Fano factor F c 12 is depicted as a function of detec- 
tion efficiencies 770 and 7712 for uq — > 00. In fact, for 
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FIG. 7: JSCF quasi-distribution P S (W , W12) of the single- 
field (Wo) and compound-field (W12) integrated intensities of 
photons for s = 0.4 (a) and the corresponding contour plot 
(b). 



Vo > ?7o.crit,mm = 0.7585 we have a sub-Poissonian be- 
havior, which does not depend on the value of 7712 (we 
suppose that 7712 = 771 = 772). This means that when 
770 exceeds its critical value, the quantum detection effi- 
ciency corresponding to the conditionally prepared state 
will not affect the nonclassical behavior of the state itself. 

Note that, from the experimental point of view, the 
nonclassical behavior of an optical state is usually ex- 
hibited by the condition R < 1. On the other hand, in 
accordance with the theoretical study presented in [28| . 
nonclassical fields satisfy the condition -Ko,i2 < 0. Ac- 
tually, the two conditions can be linked together. In 
fact, if Bo = B\2 = B, negative values of the determi- 
nant K 0,12, which testify a nonclassical character, lead 
to <[A(W - W12)] 2 ) = 2M{B 2 - |D ,i2 1 2 ) < 0. By using 
Eq. (fTT)) we thus obtain that R < 1, i.e. the sub-shot- 
noise reduction of fluctuations in the difference between 
the photons in the single and compound fields. How- 
ever, it might happen that R > 1 for Bq ^ Byi and 
^0,12 < 0: in this case the noise reduction factor cannot 
be used to check the nonclassical nature of the measured 
state. However, when this condition occurs other quan- 
tities, such as the conditional Fano factor i 7 ' c .i2, can pro- 
vide evidence of the nonclassical character of the fields. 
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FIG. 8: Conditional Fano factor -F c ,i2 as a function of the 
quantum efficiency 77 = 770 = 771 = T72 for the fixed number 
no — 3000 of the single-field photons (a) and critical quantum 
detection efficiency function of the number no of 

the single-field photons (b). If 77 > 77 cr i t then F c> r2 < 1 and 
viceversa. In graph (b), logarithmic scale is used on the x- 
axis. 



V. CONCLUSIONS 



tional photon-number distribution as well as the joint 
quasi-distribution of integrated intensities have been de- 
termined to study the nonclassical properties of the mea- 
sured fields. It has been shown that states entangled in 
photon numbers can occur in the second and third fields 
provided that a given number of photons is detected in 
the first field. In this entangled state, sum of photon 
numbers in the second and third fields equals the given 
number of photons in the first field and photon number 
in the second field (as well as in the third field) is not 
determined. The crucial role played by the quantum de- 
tection efficiencies in the production of the conditional 
state is widely discussed. 




FIG. 9: Contour plot of the conditional Fano factor -F c ,r2 as a 
function of the quantum efficiencies 770 and 7712, 7712 = 771 = 7/2 
for high values of no (no — ► 00). The values in the white area 
in the upper left corner have not been determined. F c ,i2 < 1 

for ?7o > 77o,crit,min, T)0 ,crit,min = 0.7585. 



The quantum properties of two interlinked nonlinear 
interactions generating a tripartite entangled state have 
been analyzed with special attention to the mutual cor- 
relations in the number of photons. In particular, we 
have taken into account the correlation between the pho- 
tons in the first field and the sum of photons in the other 
two. The joint photon-number distribution, its condi- 
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